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Abstract. We study an alternative to the prevailing approach to mod- 
elling qualitative spatial reasoning (QSR) problems as constraint satis- 
faction problems. In the standard approach, a relation between objects 
is a constraint whereas in the alternative approach it is a variable. 
By being declarative, the relation-variable approach greatly simplifies 
integration and implementation of QSR. To substantiate this point, we 
discuss several specific QSR algorithms from the literature which in the 
relation-variable approach reduce to the customary constraint propaga- 
tion algorithm enforcing generalised arc-consistency. 



1 Introduction 

Qualitative spatial representation and reasoning (QSR) 6} lends itself well to 
modelling by constraints. In the standard approach, a spatial object, such as a 
region, is described by a variable, and the qualitative relation between spatial 
objects, such as a topological relation between two regions, contributes a con- 
straint. For many QSR calculi, it is known that if all the constraints represent 
definite (base) relations and path-consistency (PC) holds, then this description of 
a spatial scene is consistent. If the relation is not fully specified, the correspond- 
ing constraint is a disjunction of basic constraints. By establishing PC, such a 
disjunctive constraint is refined in view of the constraints with which it shares 
a variable. A combination of PC with search over the disjunctive constraints 
decides the consistency of indefinite scene descriptions. 

We examine here an alternative constraint-based formulation of QSR. In this 
approach, a spatial object is a constant, and the relation between spatial objects 
is a variable. We call this the relation-variable approach, in contrast to the 
conventional relation- constraint approach above. Although modelling QSR with 
relation variables is not original, see |27| . it is mentioned very rarely. This fact 
surprises in view of the advantages of this approach. In particular, the following 
two important issues are tackled successfully: 

Integration. Space has several aspects that can be characterised qualitatively, 
such as size, shape, orientation. These aspects are interdependent, but no 
convenient canonical representation exists that provides a link (the role of 
time points in temporal reasoning). Spatial reasoning problems in practice 



are also not likely to occur in pure form. They may be embedded into a 
non-spatial context, or contain application-specific side constraints. 
The relation- variable approach to QSR is declarative in a strict sense and is 
thus well-suited for these integration problems. 
Systems. Typical current constraint solving platforms focus on domain reduc- 
tion, and accordingly provide convenient access to variable domains. Mod- 
ifying the constraint network, on the other hand, is usually difficult. This 
task is, however, required for enforcing PC. 

A formulation of QSR according to the relation- variable approach means that 
generic domain-reducing propagation algorithms and conventional constraint 
solving platforms can be used instead of dedicated spatial reasoning systems. 

Plan of the paper. We begin by introducing briefly the necessary constraint solv- 
ing concepts and methods, and qualitative spatial reasoning, using the example 
of the RCC-8 calculus. The next section presents in-depth the two modelling ap- 
proaches for constraint-based QSR. In the following sections, we discuss several 
aspects of space and contrast the relation- variable and relation-domain approach. 
We finally mention some new modelling options, and end with a summary. 

1.1 Constraint Satisfaction 

Recent coverage of the field can be found in |1I8I12| . 

Consider a sequence X = Xi,...,x m of pairwise different variables with 
respective domains D\, . . . ,D m . By a constraint C on X, written C(X), we 
mean a subset of D± x • • • x D rn . The arity of C is m. A constraint satisfaction 
problem (CSP) consists of a finite sequence of variables X — x%,...,x n with 
respective domains V = Di, ... ,D n , and a finite set C of constraints, each on 
a subsequence of X. We write it as (C;xi £ D\, . . . ,x n £ D n ), or shorter as 
(C; IeD). Given an element d = d\, . . . , d n of D\ X . . . X D n and a subsequence 
Y = Xi x , . . . , Xi e of X we denote by d[Y] the sequence di 1 , . . . ,di t ; in particular, 
we have d[xk] = dk- A solution to (C; X E V) is an element d £ D such that for 
each constraint C £ C on the variables Y we have d[Y] £ C. 

Constraint propagation. One method to establish satisfiability of CSPs when 
the search space is finite is systematic search for a solution. For reducing the 
search space and overall search effort, constraint propagation is often very useful; 
the principle is to replace a given CSP by another one that is equivalent with 
respect to the solutions but that is easier to solve. Constraint propagation is 
typically characterised by the resulting local consistency. The two notions most 
relevant for this paper are: 

Path Consistency (PC): A CSP of binary constraints is path-consistent [2*1| 
if for every triple of variables x,y, z 

C(x, z) = { (a, c) | b exists s.t. (a, b) £ C(x, y) and (b, c) £ C(y, z) } . 

It is assumed here that a unique constraint C(u, w) for each pair of variables 
u,w exists, and that C(u,w) — C^ 1 (w,u). 
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Generalised Arc- Consistency (GAC): A constraint C(X) is generalised 
arc-consistent [231 if f° r au x k S X and all a G D/c 

d G C{X) exists such that d[xk] = a. 

In short, every domain value must participate in a local solution. 
A CSP is generalised arc-consistent if each of its constraints is. 

For example, the CSP (x + y — z; x, y, z G {1, 2, 3}) can be reduced to 
(x + y = z; x, y G {1, 2}, z G {2, 3}) which is GAC. 

Enforcing PC means reducing constraints but not domains, whereas enforcing 
GAC means reducing domains but not constraints. 

A number of generic methods to establish GAC for a constraint are known, 
and many constraint solving systems have implementations. One example is the 
GAC-schema available in ILOG Solver [TB] , 

1.2 Qualitative Spatial Reasoning 

The topological calculus RCC-8 [2EJ is one of the best-known formalisations 
in spatial reasoning. We use it to illustrate a number of concepts. In RCC-8 
one distinguishes 8 topological relations between two regions, see Fig.^ discon- 
nected, externally connected, partially overlapping, equal, tangential proper part, 
non-tangential proper part, and inverses of the latter two. These are denoted 
DC, EC, PO, EQ, TPP, NTPP, TPPi, NTPPi, respectively; together they form a set 
that we call RCC8. 



TPP NTPP 




Fig. 1. RCC-8 relations (2D example) 



Jointly exhaustive and pairwise disjoint. Any two spatial regions are in one and 
exactly one of the RCC-8 relations to each other. 

Composition table. Considering the triple R a ,b, Rb,c, R a ,c of relations between 
regions a, b, c, one finds that not all triples of RCC-8 relations are semantically 
feasible. The consistent triples are collected in the RCC-8 composition table. It 
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contains 193 relation triples, such as (NTPP, EC, DC). Bennett |2] proved that 
compositional consistency entails global consistency: if for all triples of regions 
the relations between them respect the composition table then this topological 
scenario is consistent. 

Converse relation table. In analogy to the composition table, it is helpful to 
think of a converse relation table consisting of the 8 pairs (R, Ri) of RCC-8 
relations such that Ri is the converse of R. It contains for example (EQ, EQ) and 
(TPP, TPPi). If we agree on (EQ) for the relation of a region with itself then the 
converse relation table follows from the composition table. 

2 Approaches to Constraint-based QSR 

A spatial topological scenario consists of a set of region names denoted by 
Regions, and possibly some restrictions on the topological relation for regions 
pairs. A scenario is fully specified if for each region pair exactly one RCC-8 
relation is given. 

We examine now how scenarios can be modelled as constraint satisfaction 
problems. We continue using topology with RCC-8 as an example, but most of 
the concepts below are immediately transferable to other spatial aspects. 

2.1 Relations as Constraints 

In this conventional approach, Regions is considered to be a set of region vari- 
ables. Their infinite domain is the set of all spatial regions in the underlying 
topological space; for example, if we model 2D space then a region variable rep- 
resents a set of points in the plane. Information about the topological relation 
between two regions is expressed as a binary constraint Rel that corresponds to 
a subset of RCC8. One usually writes this in infix notation as 

constraint x Rel y where Rel C RCC8 and x,y 6 Regions. 

Such a CSP describes a possibly partially specified scenario. Whether a cor- 
responding fully specified and satisfiable scenario exists is checked by path- 
consistency and search over the relations. A PC-enforcing algorithm revises the 
constraints between regions according to the converse relation and composition 
tables of RCC-8, and search branches over disjunctive constraints. 

Establishing satisfiability of a scenario processes only the constraints, for 
compositional consistency. The variables remain unassigned. 

2.2 Relations as Variables 

Here we interpret every element of Regions as a constant. The topological rela- 
tion between two regions is a variable with a subset of RCC8 as its domain. Such 
a relation variable exists for each ordered pair of regions, and we collect all these 
variables in an array Rel. We write an individual relation as 

variable Rel[a,b] where Rel[a, b] C RCC8 and a, b £ Regions. 
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Integrity constraints. Relation converse and composition in this setting are 
captured at the constraint level. The binary constraint conv represents the con- 
verse relation table: 

conv( Rel[a, b], Rel[b 7 a] ) for all {a, b} C Regions. 
The composition table is represented by the ternary constraint comp, with 

comp( Rel[a, b], Rel[b,c], Rel[a,c\) for all {a, b, c} C Regions. 
In presence of 

Rel[a,a] — EQ for all a 6 Regions 

and a conv constraint for all pairs of different regions, one comp constraint per 
three different regions suffices. 

2.3 Comments 

By modelling the items of interest as variables and static information as con- 
straints, the relation- variable approach yields plain finite-domain CSPs in which 
the solutions (i.e., assignments) are relevant. There is a straightforward corre- 
spondence between a solution and a fully specified, consistent scenario. Obtain- 
ing the latter from a partially specified scenario amounts to the standard task 
of solving a finite-domain CSP. 

Constructing a relation-variable model means finding integrity constraints 
that embody the intended semantics. Once that has been established, the origin 
or meaning of the constraints is irrelevant. For example, a constraint solver can 
ignore whether comp represents the composition operation in a relation algebra; 
we also discuss examples below in which other restrictions on the relations must 
be satisfied. There is thus a clear distinction between specification and execution. 
The relation-variable approach is declarative in a strict sense. 

Constraint propagation. The relation- variable approach is independent of the 
particular constraint solving method. We could, however, choose a solver based 
on search and propagation, and furthermore we could choose a GAC-enforcing 
propagation algorithm. 

Path consistency in the relation-constraint approach and generalised arc- 
consistency in the relation-variable approach simulate each other. This can be 
seen by analysing, in both approaches, the removal of one topological relation 
from the disjunctive constraint a Rel b, or from the domain of the variable 
Rel[a, b], respectively. The reason in both cases must be the lack of support- 
ing relations between a, c and b, c, for some third region c; that is, compositional 
consistency. 
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Complexity. It is perhaps not surprising but useful to mention that establishing 
the respective local consistency in either approach (i.e., PC and GAC) requires 
the same computational effort. Let n denote the number of regions. Enforcing 
PC by an algorithm as the one given in |2J requires time in 0(n 3 ) [221 - F° r this, 
one assumes that one PC step, restricting a Rel c by a Rel b and b Rel c, takes 
constant time. 

Analogue reasoning entails that GAC can be enforced in constant time on a 
single comp(Rel[a, b], Rel[b, c], Rel[a, c]) constraint — observe that the three vari- 
ables have domains of size at most eight. In this way, the overall time complexity 
depends only on the number of such constraints, and is thus in 0(n 3 ). 

Previous work. Tsang |27] describes the relation-variable approach in quali- 
tative temporal reasoning, a field similar to QSR. The idea appears not to have 
caught on, however. One reason is probably that integration in temporal reason- 
ing is simpler because the canonical representation of time points on the real line 
exists. By referring to its end points, a time interval can directly be related to 
its duration or another time interval. Space, in contrast, has no such convenient 
canonical representation — but many aspects to be integrated. 

In QSR, the possibility of the relation-variable approach is mentioned occa- 
sionally in passing, but without examining its potential. For actually modelling 
and solving QSR problems using relation variables I am only aware of ^ pages 
30-33], which deals with a single aspect (topology) only. 

3 Relation Variables in Use 

An essential advantage of the relation- variable approach is that the relevant in- 
formation is available in variables. This means that linking pieces of information 
reduces to merely stating additional constraints on the variables. In that way, 
embedding a QSR problem into an application context or adding side restric- 
tions, for example, can be dealt with easily and declaratively. 

We illustrate the issue of composite models with the case of aspect integra- 
tion. 

3.1 Combining Topology and Size 

Following Gerevini and Renz J2j , we study scenarios combining topological and 
size information. We collect information about both these aspects and their link 
in one CSP. 

Let n be the number of regions. 

Topological aspect. As in Section the 

n x n array TopoRel 

of RCC-8 relation variables stores the topological relation between two regions. 
The integrity constraints convRccSi c ° m PRCC8 nee d to hold. 



6 



Size aspect. Relative size of regions is captured by one of {<, =, >}, as in [T3] . 
The 

n x n array SizeRel 

of variables stores the relative sizes of region pairs. The converse relation and 
composition tables are straightforward; the integrity constraints are 

conv Sizc = { (<, >), (=,=), (>,<)}, and 
comp Sizc = { (<, <, <), (<, =, <), . . .} (13 triples). 

Linking the aspects. The topological relation between two regions is depen- 
dent on their relative size. A table with this information is given in it 
contains rules such as the following: 

TopoRel[x, y] = TPP implies SizeRel[x,y] = (<), 

SizeRel[x,y] = (=) implies TopoRel[x 7 y] G {DC, EC, PO, EC?}. 

In |13j . these rules represent a meta constraint. Here, we infer the linking con- 
straint 

linkxopo&Sizo - { (TPP, <), (DC, =), . . .} (14 pairs) 
which is to be stated as 

linl<Topo&Sizc( TopoRel[a,b], SizeRel[a,b]) 
for all regions a, b. 

Example. Let us pick up the combined scenario from |13l p. 14]. Five regions, 
denoted by {0, . . . , 4}, are constrained by 

TopoRel[Q, 2] G {TPP, EQ} SizeRel[0, 2] G {<} 

TopoRel[l, 0] G {TPP, EQ, PO} SizeRel[3, 1] £ {<, =} 

TopoRel[l,2] G {TPP, EQ} SizeRel%A] e {<,=} 
TopoRel[A,i] £ {TPP, EQ} 

Independently, the topological and the size scenarios are consistent while the 
combined scenario is not. It is pointed out in that naive propagation schedul- 
ing schemes do not suffice to detect inconsistency. 

A formulation of this scenario as a combined topological & size CSP in the 
relation-variable approach is straightforward. The resulting CSP can be entered 
into a constraint programming platform such as ECL l PS e [2H|- ECL^PS 6 is 
focused on search and domain-reducing propagation; in particular, it offers a 
GAC-enforcing propagation algorithm for user-defined constraints. Given our 
CSP in ECL l PS e , solely executing GAC-propagation for all constraints yields 
failure, which proves that this CSP is inconsistent. □ 
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For the same purpose but within the relation-constraint approach, Gerevini 
and Renz proposed a new algorithm called Bipath-CONSISTENCY ^3] ■ Its princi- 
ple is the computation of path-consistency for both types of relations in an inter- 
leaved fashion while taking into account the interdependency. The linkxopo&sizo 
constraint is in essence treated as a meta constraint on the algorithm level. More- 
over, the Bipath-CONSISTENCY algorithm fixes in part the order of propagation. 

The relation-variable method, on the other hand, is declarative; all informa- 
tion is in the five types of constraints. They are handled by repeated, interleaved 
calls to the same GAC-enforcing algorithm. The actual propagation order is ir- 
relevant for the result. 

Bipath-CONSISTENCY is restricted to combining two types of relations (e.g., 
two aspects of space). In contrast, the relation- variable approach is composi- 
tional in the sense that adding a third aspect, such as morphology 7 or orien- 
tation, is straightforward. It amounts to formulating integrity constraints (e.g., 
conv, comp), linking constraints to each of the already present aspects, and a 
constraint linking all three aspects. Some of these constraints may be logically 
redundant. 

3.2 Combining Cardinal Directions and Topology 

In orientation, another important aspect of space, one studies the relation of 
two objects, the primary and the reference object, with respect to a frame of 
reference. It is thus inherently a ternary relation, but by agreeing on the frame 
of reference, a binary relation is obtained. 

The binary relation approach is realised in the cardinal direction model [§], 
based on the geographic (compass) directions. Points as well as regions have 
been studied as the objects to be oriented. The point-based models can be cast 
in the relation- variable approach analogously to topology, Section 1^1 For in- 
stance, Frank distinguishes the jointly exhaustive and pairwise disjoint re- 
lations N, NW, W, . . . for points; denoting North, Northwest, West, and so on. 
Ligozat [201 gives a composition table. 

Orienting regions. Goyal and Egenhofer and Skiadopoulos and 

Koubarakis [2U study a more expressive model, in which the oriented objects 
are regions. The exact shape of the primary region is taken into account, and 
a ninth atomic relation B exists, describing overlap of the primary region and 
the axes-parallel minimum bounding box of the reference region. Sets of the 
atomic relations are then used to describe directional information. In this way, 
for example, the position of South America for an observer located in Ecuador 
can be fully described by the set {B, N, NE, E, SE, S}. In contrast, the position of 
Ecuador with respect to South America is just {B}. 

Relation variables for directional information are thus naturally set variables: 
they take their value from a set of sets of constants, unlike relation variables for 
topology and size whose domain is a set of atomic constants. 
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For each pair a, b of regions, the direction is a relation variable 



DirRel[a, b] e P(Dir) where Dir = {B, N, NW, . . . , NE}. 



V denotes the power set function. 



Integrity constraints. A restriction on the set values that DirRel[a, b] can 
take arises if a, b are internally connected regions, which is often assumed. Only 
218 of the 512 subsets of Dir are then semantically possible. This knowledge 
can be represented in a unary integrity constraint, which for example allows 
{N, NE, E} but excludes {N,S}. The usual integrity constraints comp and conv 
can be derived from studies of composition ^Zjg\ and converse jS] (but it is outside 
of our focus whether these are the only integrity constraints needed). 



Integration with topology. Let us briefly consider linking directional infor- 
mation to topology. The relevant knowledge could be expressed by rules as 

TopoRel[x,y] <E {EQ, NTPP, TPP} implies DirRel[x,y] = {B}, 
TopoRel[x, y] £ {NTPPi,TPPi} implies DirRel[x, y\ D {B}, 

from which a constraint linkxopo&Dir can be defined. It is to be stated as 

linkxopo&Dir( TopoRel[a, b], DirRel[a, b] ) 

for all regions a, b. We now have some components of a combined cardinal direc- 
tions & topology model. It can be given to any sufficiently expressive constraint 
solver, which in particular would provide constraints on set variables. 

Constraint solving with set variables is discussed in |14j . Many contemporary 
constraint programming systems support set variables. 



3.3 Cyclic Ordering of Orientations with Relation Variables 

From the several formalisations of orientation information with an explicit frame 
of reference, let us examine the approach of Isli and Cohn to cyclic ordering of 2D 
orientations [T2|. Here, the spatial objects are orientations, i.e. directed lines. At 
the root of the framework is the qualitative classification of the angle a = <(a, b) 
between the two orientations a and b by 



Or(or) = 



(equal) if a = 0, 

(left) if < a< 7T, 

o (opposite) if a = n, 

(right) if 7r < a < lix 



into the jointly exhaustive and pairwise disjoint relations e, l,o, r. See Fig. |2]for 
an illustration. For three orientations a, b, c, we now consider the triple 
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Fig. 2. The relations e, I, o, r of a pair of orientations 






Fig. 3. The Cyc relations Irl, lei, and rol of a triple of orientations 



(Or(<(M)), 0r(<(c,6)), Or(<(ca)) >. 

Of all 4 3 triples over {e, I, o, r}, only 24 combinations are geometrically possible. 
We denote this set by Cyc. Fig. [3| shows three of its elements. 

Such cyclic ordering information can be expressed within the relation variable 
approach in an array CycRel, which in particular is ternary. We have thus a 
relation variable 

CycRel[a, b, c] G Cyc with Cyc = {Irl, orl, . . . , rle} 

for every three orientations a, b, c. The integrity constraints here are 

conv( CycRel[a,b, c], CycRel[a 1 c 1 b\) 1 

comp( CycRel[a, b, c], CycRel[a, c, d\, CycRel[a,b,d}), 

and a new constraint 

rotate( CycRel[a, b, c], CycRel[c, a, b] ). 

Details and definitions can be found in |19j . 

Working within the relation constraint approach, Isli and Cohn construct a 
new algorithm called s4c that enforces 4-consistency on the ternary relation 
constraints that correspond to CycRel. They are able to prove that this algorithm 
decides consistency, i.e., 2D geometric feasibility, of fully specified scenarios. The 
s4c algorithm uses exactly the information that we represent in the conv, comp 
and rotation constraints. Consequently, we can conclude that in our relation 
variable model these constraints guarantee geometric consistency. 
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We hypothesise further that s^c in the relation constraint model propagates 
at most as much information as a GAC-enforcing algorithm does in our relation 
variable model. Intuitively, this should be clear: every possible reduction of a 
disjunctive constraint in the relation constraint model corresponds to a domain 
reduction of a relation variable in our model. 

3.4 Combining Cardinal Direction with Relative Orientation 

Isli |17I18| studies the problem of exchanging information between a cardinal 
direction model for pairs of points as in Section 13.21 and a relative orientation 
model for triples of points, derived from Freksa and Zimmermann's formalisa- 
tion ^D]. This problem is again similar to combining topology and size, Sec- 
tion Isli works with the relation-constraints and proposes a new algorithm 
for this integration issue. 

We formulate a relation-variable model. The cardinal direction subproblem 
can straightforwardly be expressed in this approach; we omit the obvious details 
here. The relative orientation subproblem leads to a model similar to that of 
orientations in the preceding section; in particular, it is based on a ternary 
array. The arrays in the combined model are: 

n x n x n array ROrientRel, and 
n x n array CDirRel, 

if we assume n points. 

For linking the two models, Isli JS| devises functions for both directions of 
the information transfer. They can be transformed into the two constraints 

linkcD->R.o( CDirRel[a, b], CDirRel[b, c], ROrientRel[a, b, c] ), 

linkcDi— ro( ROrientRel[a, b, c], C'DirRel[a,b], CDirRel[b, c], CDirRel[a, c] ). 

For the relation-constraint model it is necessary to treat the information in 
linkcD— i-ROj li n kcD<-RO a s meta-constraints, embedded inside an algorithm that 
moreover integrates sJ^c of |19j and a path-consistency algorithm. 

Using relation variables, it suffices to state the constraints and provide a 
generic GAC-enforcing algorithm. Also, for a given triple of points, the first 
constraint linkcD->RO should just be the restriction of the second constraint 
linkcD<-RO in which the variable CDirRel[a,c] is projected away. The former 
constraint is then redundant, and we just need one constraint 

linkcD&Ro( ROrientRel[a, b, c], CDirRel[a,b], CDirRel[b, c], CDirRel[a, c] ). 

On the grounds that both the relation-variable and the relation-constraint 
approach are based on the same semantic information, for one embedded in an 
algorithm, for the other in constraints, we conclude that both accept exactly the 
same point configuration scenarios. 
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4 Extensions 

Variables ranging over spatial objects. In the relation-variable model, spa- 
tial objects are denoted by constants. An object variable, whose domain is the set 
of object constants, has thus a different meaning than in the relation-constraint 
approach. This issue is best demonstrated by an example. Suppose we wish to 
identify two regions among all given regions such that 

— the first is smaller than the second, and 

— they are disconnected or externally connected. 

We use topological and size information, formalised as in Section Em so we have 
arrays SizeRel and TopoRel recording the qualitative relations. Let Regions be 
the set of the n region constants. We define the 

region variables x\ , X2 

whose domain is the set Regions, and constrain them by 



C\ is a constraint on the variables X\, xi and on all size relation variables in the 
array SizeRel. Namely, region constants T\,t2 € Regions must be assigned to 
xi,X2 such that the size relation variable SizeRel[r\, r 2 ] is assigned a '<'. 

We call such constraints, in which arrays are indexed by variables instead 
of constants, array constraints. They are a generalisation of the better-known 
element constraint, which corresponds to a one-dimensional array indexed by 
a variable. Constraint propagation to establish GAC for array constraints is 
studied in fQ. The constraint programming system ILOG Solver |lfi| accepts 
and propagates array constraints. 

Reasoning about spatial change. It is not difficult to augment a relation- 
variable model with temporal information. It suffices to add a new time index to 
each array of qualitative relations, and to link the new time-annotated scenarios 
appropriately. We extend Rel from a binary to a ternary array such that 



is a variable specifying the relation between the spatial objects a and b at time 
t. Suppose we view time as linear and discrete, such that only atomic relational 
changes can occur between subsequent time points. We can specify these atomic 
changes (the so-called conceptual neighbourhood) by pairs of qualitative rela- 
tions and define accordingly a new binary constraint neighbour. For example, 
the pair (DC, EC) in the constraint neighbour Topo indicates that the topological 
relation disconnected between two regions may change in one time step to ex- 
ternally connected. The neighbour constraint is then stated on all variable pairs 
(Rel[a, b, t],Rel[a, b, t']) where t directly precedes t' temporally. 



SizeRel[x\,X2] = (<), 
TopoRel[xi, xa] e {DC, EC}. 



(Ci) 
(C 2 ) 



Rel[a,b,t] 
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5 Summary 



We have presented an alternative formulation of qualitative spatial reasoning 
problems as constraint satisfaction problems. Contrary to the conventional ap- 
proach, we model qualitative relations as variables. Uncertain relational infor- 
mation is naturally expressed by variables with domains; consistency of this in- 
formation is naturally expressed by static constraints. The propagation of these 
constraints is a well-understood issue in research on constraint programming, 
and corresponding generic algorithms are provided by many constraint solving 
systems. 

While the principle of the relation-variable approach is not new, the advan- 
tages of applying it to QSR, especially for integration tasks, have so far very 
rarely been realised. We have argued that several algorithms that are custom- 
designed for integrating spatial aspects become unnecessary if a relation- variable 
model and a generic GAC-establishing constraint propagation algorithm is used: 
the Bipath-CONSISTENCY algorithm of ! :'> , the sJ^c algorithm of ^5], the algo- 
rithm combining s4c and a path-consistency algorithm of JSj. We have shown 
how the relation-variable approach can accommodate composite qualitative re- 
lations as investigated in 5 26. with the help of set variables and constraints. 
We have indicated that extending or combining a relation-variable model often 
consists mainly in defining appropriate constraints, contrary to what is the case 
in the relation-constraint approach where new algorithms must be designed. 

Finally, we remark that the strictly declarative model that results from using 
relation-variables can be solved by any sufficiently expressive solver of CSPs. 
This includes typical CP systems based on search and propagation, but also for 
example solvers based on local search. 
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